Charge oscillations in Quantum Dots: Renormalization group and Hartree method 

calculations 



in 
o 
o 

(N 

> 
O 

00 



c3 



i 

C 

o 
o 



> 

a\ 
o 
oo 
o 
^t- 
o 

-I— > 



i 

C 

O 

o 



Michael Sindel, 1 Alessandro Silva, 2,3 Yuval Oreg, 3 and Jan von Delft ^ 

1 Physics Department and Center for NanoScience, LMU Munchen, 80333 Miinchen, Germany 
2 Center for Materials Theory, Department of Physics and Astronomy, Rutgers University, Piscataway, NJ 08854, USA 
3 Department of Condensed Matter Physics, The Weizmann Institute of Science, Rehovot 76100, Israel 

(Dated: February 2, 2008) 

We analyze the local level occupation of a spinless, interacting two-level quantum dot coupled to 
two leads by means of Wilson's numerical renormalization group method. A gate voltage sweep, 
causing a rearrangement of the charge such that the system's energy is minimized, leads to oscil- 
lations, and sometimes even inversions, in the level occupations. We find that these oscillations, 
qualitatively understandable by a simple Hartree analysis, are generic and occur in a wide range of 
system parameters. By allowing a relative sign in one tunneling matrix element between dot and 
leads, we extend our findings to more generic models. Experimental applications and the qualitative 
effect of spin are discussed. 

PACS numbers: 73.63.Kv, 73.23.Hk, 72.10.-d 



I. INTRODUCTION 

The Coulomb blockade (CB)i is one of the most ba- 
sic and fundamental phenomena in quantum dot (QD) 
physics: according to the standard single-electron tun- 
neling picture for transport through small QD's, elec- 
trons can only pass through the dot one by one if the 
Coulomb charging energy U is the dominant energy scale. 
Likewise, if the gate voltage (applied to a plunger gate 
near the QD) is changed, resulting in a rigid shift of the 
single-particle spectrum of the QD relative to the Fermi 
energy of the leads, the single-particle levels of the QD 
are naively expected to be filled "one by one" . Many ba- 
sic single-electron devices and QD-based qubit proposals- 
are based on this simple picture. 

In this paper we point out that under quite generic 
conditions, the charging process is actually more compli- 
cated, and the occupation functions for the single-particle 
levels can show some rather complex, non-monotonic be- 
havior as a function of gate voltage, deviating consider- 
ably from the standard CB picture ["one by one" filling]. 
This complex behavior turns out to be rather generic for 
a QD coupled to leads that can supply the QD with elec- 
trons, and originates from the competition between the 
QD-lcads coupling T, and the intrinsic energy scales of 
the QD, namely its charging energy U and levels spacing 
A. 

In order to study this competition, it suffices to con- 
sider a very simple model: a spinless interacting two-level 
Anderson model (2LAM), consisting of a lower and an 
upper QD level (ee, e u ) with level spacing 3 . A = e u — eg. 
Analyzing this model using both the numerical renor- 
malization group (NRG)ii and a self-consistent Hartree 
approximation, we study in detail the evolution of the 
occupation of the single particle levels as a function of 
gate voltage at T — for various values of the QD pa- 
rameters (r, A, U). The generic picture of the charging 
process emerging from this analysis is the following: for 
any finite coupling to the leads T (-^C A), sweeping the 



levels towards the Fermi level of the leads (by tuning the 
gate voltage) causes the occupations of both the lower 
and the upper level to increase at comparable rates. This 
process continues until one level takes over and becomes 
more occupied than the other. At this point the electron 
that occupies this level electrostatically repels the other 
level, pushing up its energy and thereby emptying it. As 
a result the occupation of the other level performs an 
oscillation as the gate voltage is swept. The naive QD 
charging scheme, in which every step of the CB stair- 
case is associated with the filling of only one single par- 
ticle level in the QD, is only achieved when these charge 
oscillations are small, i.e. for r <C A. Below, we dis- 
cuss in detail the physics of "charge oscillations" and the 
dependence of their amplitude and form on the system 
parameters. In particular, we discuss under which condi- 
tions the amplitude of these oscillations can be made so 
large as to cause a significant "occupation inversion" , i.e. 
a situation for which the occupation of the lower level is 
smaller than that of the upper level. It is important to 
mention that more general models, e.g. including spin, 
show effects similar to those described above. 



The rest of the paper is organized as follows. In Sec- 
tion II we describe in detail the model we consider. In 
Section III, we present the result of the NRG calculation 
for the evolution of the occupation of the two levels as 
a function of gate voltage [Section III A], and a simple 
and qualitative understanding of the emerging picture in 
terms of an Hartree approach [Section III B]. In Section 
IV, we present a detailed analysis of the phenomenon 
of occupation inversions, by studying the 2LAM in the 
case where the two levels are unequally coupled to the 
leads. Finally, in Section V we summarize our conclu- 
sions, and discuss the robustness of the effects obtained 
when including spin, and possible experimental tests for 
our predictions. 
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FIG. 1: Schematic depiction of the model. The sign s between 
the tunneling amplitudes VgR and Vil determines whether the 
two dot levels couple to the same (symmetric) channel (s = 
+1, upper mapping), or to two different channels (s = — 1, 
lower mapping), with strength Vi = V2Vi, i = l,u. 



namely n u (e + e*) = 1 — ng(e* — e), independent of A and 
r. If 7^ r„, this symmetry is broken for both s = ±1. 



III. CHARGE OSCILLATIONS 

In this section, we study the evolution of rii (i = 1, 2) as 
a function of the gate voltage e, and discuss the physical 
origin of non-monotonic occupation of the two levels. We 
shall first present the results of our NRG 1 - 1 calculations 
for m, and then use a simple Hartree analysis to gain 
some insight into the NRG results. 

A. NRG calculations 



II. MODEL 

We consider a spinless 2LAM with Hamiltonian H = 
Hd+Ti-i +Hid- (The spinful case will be considered briefly 
in Section^]) Two leads, identical, noninteracting and in 
equilibrium, are described by Hi =^2 ka ^k a c\ a Ck a , where 
c ka crea tes an electron with energy in lead a — L, R. 
The isolated QD is described by 

U d = } J e^X + Unth u , (1) 

i— u,£ 

where d\ creates an electron in the QD in level e; 
(i = £,u), measured w.r.t. the Fermi level defined by 
the leads, hi = d\d i is the number operator, and U 
is the charging energy, which we fix at U = 0.2D 
throughout this paper, 2D being the bandwidth. Fi- 
nally, the tunneling between the QD and the leads is de- 
scribed by Hu = Y, kia {Vkia, c ka d i +h.c). We consider k- 
independent tunneling matrix elements Vki a = Via , which 
are L-i?-symmetric in magnitude 4 , V u l = V u r = V u and 
ViL = sVgR = Vg, but with a possible relative phase 
s — exp(i0) between the L and R matrix elements of 
the lower level^ffiiS. Time reversal symmetry implies 
4> = 0, 7T (hence s — ±1). The corresponding bare level 
widths are = 2irpVj 2 , where p is the density of states 
in the leads. 

The two possible choices for s — ±1 lead to two dis- 
tinct models (see Fig.QJ: (i) for s — +1, both local levels 
couple to the same channel, namely the symmetric linear 
combination of the left and the right lead (ckL + CkR,), 
while the antisymmetric combination (ckL ~ CkR.) decou- 
ples completely; (ii) for s = —1, the upper and lower 
local levels couple to different channels, namely to the 
symmetric or antisymmetric combinations, respectivelyS,. 



We start our analysis by considering equally coupled 
levels, = r n = r, and use the NRG 11 to calculate 
the e-dependence of r^. Naively one may expect that 
if the QD is initially empty, the QD levels get occupied 
monotonically one by one as e is decreased, the usual 
CB behavior. In other words, first the occupation of the 
lower level would be expected to increase monotonically 
as eg crosses the Fermi level, and subsequently n u would 
increase as e u + U approaches it. However, our NRG re- 
sults [Fig. |2] show that this intuitive picture is valid only 
if the coupling to the leads is much smaller than the dot 
level spacing, T <C A. In particular, when r > A, rig 
and n u show a non-monotonic e-dependence, character- 
ized by charge oscillations of n u (or ng) when the lower 
(or upper) level crosses the Fermi leveti^. The oscillation 
in n u occurs because as soon as the lower level begins to 
be occupied significantly, the system can gain charging 
energy by additionally filling the lower level and emp- 
tying the upper level (an analogous argument works for 
ng). 




We shall denote the ground state expectation value of 
the occupation of level i by ni — {hi), and parameter- 
ize the gate voltage by the average bare level position 
e = (eg + e u ) /2. For Tg — T u = T, this parameteriza- 
tion reveals particle-hole symmetryifi around e* = —U/2, 



FIG. 2: NRG results for the occupation of the lower [ng, 
(curves 1' - 3')] and the upper level [n u , (curves 1 - 3)] for 
fixed A and different values of F, for (a): s = +1 and (b): 
s = —1. The dotted lines indicate where the lower/upper 
level crosses the Fermi level, at e = A/2 and e = —U — A/2. 
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To explore how strongly these charge oscillations vary 
with r, A and how they are affected by the sign of s, we 
show in Fig. the behavior of rij(e) for variable T and 
s = ±1, keeping A fixed (at 0.2U). In the limit T/A < 1 
level £ becomes occupied rather suddenly (curve 1') when 
it crosses the Fermi level at e = A/2, and similarly for the 
upper level at e = — U — A/2 (curve 1). In addition to this 
typical CB behavior, we observe, even for the smallest 
r considered (r = 0.1A), a tiny non-monotonicity or 
charge oscillation in n u (ng) roughly at that e where the 
occupation of the lower (upper) level increases sharply 
from to 1. A gradual increase of T towards r/A w 1 
results in a strengthening of these charge oscillations. In 
the limit T/A ^> 1, Fig. E2c) and (d), the monotonic 
increase in the occupation is recovered, though the actual 
dependence of rij on e depends strongly on the sign of s. 



B. Hartree approach 

A simple and qualitative understanding of the charge 
oscillations observed in our NRG results can be obtained 
in the framework of a self-consistent Hartree approxima- 
tion (scHA), by studying the evolution of the Hartree 
levels as function of e. This scheme accounts for the 
interaction by replacing the bare levels of the nonin- 
teracting problem, with the corresponding Hartree levels 



Uut, 



(2) 



where i = i/u if i = u/i. By integrating out the leads, 
we obtain the effective noninteracting dot Hamiltonian 



H 



err 



-WT U TA 5 s 



(~iy/T u Ti) <5 s .+i 
e¥ - iTg 



(3) 



The corresponding retarded dot Green's function, defined 
as Q${t) = -£0(t)({di(t),dt(O)}}, can be obtained ex- 
actly for both values of s, by solving the matrix equation 
Q (u>) = (u> — Ti-cfi) 1 ■ To finally obtain the Hartree ap- 
proximation for the interacting Green's function, one has 
to self-consistently calculate the average level occupation 
Tii(e), using the T = relation 



did lra,G u (w, e) 



(4) 



Since the self-consistent Hartree equation Q may have 
more than one solution, a criterion is needed to pick the 
correct one. To this end, we note that, for given e, the 
system adjusts its local level occupations n u and ng such 
that its total free energy F s (n Ul ng) is minimized. Within 
the scHA approach, F s can be obtained by integrating 
Eq. (0J, so that the conditions for F s to be extremal, 
dF s /dni = 0, reproduce Eq. we should then pick 
that solution of Eq. for which the extremum is a global 
minimum of F s . In the case of two nearly degenerate 
minima, the scHA neglects the possibility of tunneling 



between them, and a different approach has to be consid- 
ered. Nevertheless we find, somewhat unexpectedly, that 
in the case of exactly degenerate minima, e.g. A = 0, 
an average over the minima reproduces the NRG results 
rather accurately [see T/A = oo curve in Fig. E£d)]. 

For simplicity we implemented this strategy explicitly 
for s = — 1, but not for s = +1, since for the latter TCf^ is 
not diagonal, which makes the determination of F+\ very 
tedious. For s = — 1, Eq. Q yields the condition n^(e) = 
^ — — arctan{(ej + Un^) /Fj} , and the corresponding free 
energy has the form 



F_x = Un e n u + 



tii^i log (sin7mi 

7T 



(5) 



Fig-|2|compare NRG with corresponding scHA results for 
7ii(e). For s = —1, we minimized [Eq. Ip)l] and find 
remarkably good agreement between NRG and scHA. For 
s = +1, for which we did not determine -F+i, we show the 
results of a "naive scHA", obtained by simply plotting 
a numerical solution of Eq. Q and "hoping" (without 
checking) that it is the correct one. Clearly, the results 
so obtained cannot be trusted on their own merit; we 
present the naive scHA results nevertheless, to illustrate 
precisely this point: indeed, in Fig. |^a,c) [for T/A > 1] 
they do not agree well with NRG results. 

Fig. |3Jc) and (d) include a special situation, namely 
7=1 and A — 0, for which both eg = e u and ef = e^ t . 
This causes a sudden jump for s = +1 in n u , but none 
for s = -1 [cf. dashed lines in Fig. |3{c) and|^d), re- 
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FIG. 3: Comparison of NRG and scHA results for n u (e), for 
fixed A and variable F [(a), (b)] or fixed V and variable A 
[(c), (d)]. The naive scHA used for s = +1 works well for 
T/A < 1. Insets in (b) and (d): Hartree levels ef (circles) 
and (stars), calculated via the scHA from Eq. JJ}> for 
r/A = 0.1 in (b) and 0.4 in (d). For T/A = 00 the Hartree 
levels are degenerate (not shown). The arrow in inset (d) 
indicates the value of eo, the local minimum of t u . 
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spectively] . To understand why, note that for A = and 
s = +1 the odd local combination (d u — dg) decouples 
from the leads; thus, its width is zero and hence its oc- 
cupation increases abruptly when its energy drops below 
zero. On the other hand, for s = — 1 the occupation in- 
creases gradually, since the width of the odd combination 
is comparable to that of the even one (d u + dg). A similar 
argument explains why for small but non-zero A/r (odd 
level almost decoupled for s — +1) as in Fig. |2 curves 
3 and 3', the charge oscillations are still observable for 
s = +1 but not for s = — 1. 



IV. UNEQUAL COUPLING AMPLITUDES 
LEADING TO OCCUPATION INVERSION 

In addition to providing a simple way to compute phys- 
ical quantities, the scHA, and in particular the concept of 
Hartree levels [see insets of Fig. EHb,d)], may be used to 
qualitatively understand how the physics of the charge 
oscillations depends on the various system parameters. 
Suppose that both Hartree levels are swept downwards, 
starting from e well above the Fermi level. When the 
lower level comes within Tg of the Fermi level, it begins 
to fill up and the upper Hartree level is pushed up by 
U , causing a charge oscillation in n u . The latter will be 
stronger the larger n u was before the oscillation, i.e., the 
larger the width (T u ) of the upper level, and the lower the 
value e^(eo) of the upper Hartree level at its local mini- 
mum, say eo (cf. Fig. © . Indeed, if e^(e ) < r„, then the 
upper level achieves a rather significant occupation be- 
fore full occupation of the lower level (and corresponding 
emptying of the upper one), implying an increase in the 
amplitude of the corresponding charge oscillation. More- 
over, since e^(eo) is also the lower the more suddenly the 
lower level gets filled, a smaller 1^ also strengthens the 
charge oscillations. Thus, strong charge oscillations can 
be obtained quite generally by allowing 7 = r„/r^ 7^ 1. 
The above argument implies that n„-oscillations are en- 
hanced for 7 3> 1; by an analogous argument, with I <-> u, 
rig- oscillations are strengthened for 7 <C 1. 

It appears that the oscillations are so strong that when 
7^1, the lower and the upper Hartree levels might 
actually cross each other (see inset Fig. |5J| , leading to an 
inverted occupation [Figs. 14151 and IS], 

Since the bare energy levels are separated by the level 
spacing A, the conditions Un > A and max (Tg, T u ) > A 
must be met to achieve such an occupation inversion. 
Figs. 0] to show how the asymmetry (7 7^ 1) of the 
couplings affects the occupation of level I and u (dashed 
and solid lines) both for s = ±1, leading to an inversion 
of the occupation within a certain range of e. 

Our discussion of occupation inversion generalizes a 
recent related study by Silvestrov and Imryii In an at- 
tempt to understand the origin of repeated and abrupt 
phase lapses observed in the transmission phase of a mul- 
tilevel QDii, they studied a multilevel model consisting 
of one level (say u with coupling T u ) strongly coupled 




FIG. 4: Allowing 7 = Y u /Yi to be 7^ 1 results in an inversion 
of the occupation for a certain range of e. The scHA results 
(n„: crosses, n^. boxes) for s — —1 (where F—i is known) 
agree well with NRG results (n u : solid, nf. dashed lines). For 
7 < 1 the Hartree levels cross near the left CB peak, implying 
an occupation inversion below the corresponding crossings. 



S=+1 S=-1 




FIG. 5: Same as Fig. El but for 7 > 1. Now the Hartree levels 
cross near the right CB peak (see inset (b) for 7 = 4), implying 
an occupation inversion above the corresponding crossings. 
Notice the differences in the shape of the occupation curves 
between the s — 1 case and s = — 1 case. 



to the leads and at least one additional weakly coupled 
level (say I with coupling I^). Ref. 0(a) considered 
the limiting case r„ > A > Y t — > (Ref. [13(b) also 
studied finite but small Te and spin), and compared the 
energies of the configuration (n u ,ni) = (1,0) to that 
of {n u ,ni) = (0,1) in second order perturbation the- 
ory in the tunneling. Their results indicate that the 
system is able to sustain an occupation inversion until 
e u = -U/ [exp (27rA/r„) + 1] w -C7/2. Although we ap- 
proach this problem from a different angle, i.e., we either 
solve it exactly by NRG or first solve the tunnel-coupling 
exactly and then treat the interaction self-consistently, 
the inversion range found in Ref. coincides with the 
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FIG. 6: (a) NRG results for the occupations with s = 1 
and 7 = 30. (b) A comparison between the scHA results 
(n» crosses, nt boxes) and for the NRG results (n„ solid, ng 
dashed lines), for s = —1 (where F-\ is known) and 7 = 30. 
The agreement between the NRG and the Hartree approxi- 
mation is remarkably good and does not depend on 7 

results 1 - of this article. 

Our analysis indicates that the example of Ref . [l3| for 
occupation inversion is a very special case of a much more 
general phenomenon whose strength depends on 7: as 7 
is increased from 1 (where no occupation inversion oc- 
curs), (i) the range of gate voltages in which inversion 
occurs increases, with the inversion point moving towards 
the middle of the CB valley; and (ii) the maximal value 
reached by n u right before the inversion increases gradu- 
ally towards 1, i.e., the effect becomes more pronounced. 

V. SUMMARY 

In this paper, we have studied the gate voltage depen- 
dence of the occupation of a spinless two-level Anderson 
model for the generic case of a relative sign s in the tun- 
neling amplitude. We found a non-monotonic behavior in 
the occupation of the local levels, due to charging effects 
between electrons within the QD, and explained this ef- 
fect in the framework of a self-consistent Hartree approx- 
imation. Remarkably, the occupations of the upper and 
lower levels can even be inverted if the level-to-lead cou- 
plings are sufficiently asymmetric. Even though we fo- 
cused onT = throughout this paper we expect the cal- 
culated behavior to persist as long as T < min{r, A, U}. 

The inclusion of spin in the 2LAM, though making the 
problem more complex, does not change the qualitative 
results presented above in a large region of parameter 
space. To illustrate this, we show in Fig. [7| the total oc- 
cupation of the upper and lower levels obtained by NRG 
for the spinful model. The level crossing persists and ad- 
ditional oscillations are observed due to the possibility to 



put two electrons of opposite spin within each level. 




E/U 

FIG. 7: NRG results in the presence of spin (with the rela- 
tive sign s = —1). The parameters are identical to those in 
Fig. The occupation crossing persists in the spinful case. 
The additional oscillations in the occupations are observed 
due to the possibility to put two electrons of opposite spin 
with in each level. 



Our predictions of the non-monotonicity of the charg- 
ing of a 2LAM are experimentally relevant (in spirit, if 
not in detail) for any quantum dot system containing 
orbital levels that are "nearly degenerate", in the sense 
that their spacing is smaller than the level widths. One 
way of realizing the specific models studied here would be 
to use two capacitively coupled quantum dots, each with 
large level spacings, associating their topmost not-fully- 
occupied levels with e u and eg, and using a large magnetic 
field to lift the Zeeman degeneracy of each. The way in 
which the charges on these dots evolve with gate voltage, 
i.e., the evolution of rii(e), could then be measured ex- 
perimentally using QPC's serving as extremely sensitive 
charge sensors, see e.g. Fig. 1 (a) of Ref. [TlJ In addi- 
tion, we expect that the occupation of the levels will be 
reflected in other properties of the system such as the 
transmission phased We leave that, however, for future 
studies. 
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